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Abstract. The IC 348 young star cluster contains more than 300 confirmed members. It is embedded in the Perseus molecular 
cloud, making any clustering analysis subject to an extinction bias. In this work, we derive the extinction map of the cloud and 
revisit the content of IC 348 through a statistical approach that uses the 2MASS data. Our goal was to address the question of 
the completeness of IC 348 and of young clusters in general. We performed a combined analysis of the star color and density 
in this region, in order to establish the surface density map of the cluster. We reached the conclusion that IC 348 has structures 
up to 25' from the cluster center, and we estimate that about 40 members brighter than K s = 13 mag are still unidentified. 
Although we cannot use our statistical method to identify these new members individually, the surface density map gives a 
strong indication of their actual location. They are distributed in the outer regions of the cluster, where very few dedicated 
observations have been made so far, which is probably why they escaped previous identification. In addition, we propose the 
existence of a new embedded cluster associated to the infrared source MSX6C G160. 2784-18.4216, about 38' south of IC 348. 
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1. Introduction 

Young stellar clusters concentrate stars in their early stages, 
making them the privileged location for studying the initial 
conditions of star formation from the molecular cloud frag- 
mentation to the evolution on the main sequence. Such a large 
and ambitious enterprise requires the acquisition of multi- 
wavelength data from the millimeter to the X-ray spectral 
range, in order to probe the different physical phenomena that 
occur during star formation. The identification of cluster mem- 
bers leads to constrain the cluster history, the initial mass func- 
tion, and the dynamical evolution. Although the ultimate goal 
is to take a complete census of the clusters where each object 
is actually given a membership status, the individual identifica- 
tion is not easy to obtain and is not always necessary. A statis- 
tical identification could also lead to information on the cluster 
morphology and the initial mass function. 

The young cluster IC 348 has been investigated regularly 
for decades and is still a very popular target. Its brightest mem- 
ber is the B5 star BD +31 643, which is part of the Per OB2 
association (Blaauw, 1952). A distance of 320 pc is usually 
adopted for the cluster although it is still controversial (Cernis, 
1993; Herbig, 1998; de Zeeuw et al., 1999). Herbig (1954) 
identified 16 members characterized by their Ha emission, and 
the current number is about 300, with a mass and age distribu- 



Send offprint requests to: L Cambresy 
e-mail: cambresy@astro . u-strasbg . fr 



tion peak at 0.2 M Q and 2 Myr, respectively (Luhman, 1999; 
Carpenter, 2002; Luhman et al., 2003; Muench et al., 2003; 
Preibisch et al., 2003; Preibisch & Zinnecker, 2004; Luhman 
et al., 2005b). Muench et al. (2003) showed that the extinc- 
tion reaches values as high as 20 mag toward the cluster line 
of sight, the cluster itself being only partially embedded in the 
Perseus molecular cloud with an average extinction of 4.5 mag 
(Lada & Lada, 1995). The extinction will introduce a bias in 
the apparent member distribution, as long as the cluster census 
is incomplete. It is a general characteristic of young clusters to 
still be associated with their parent molecular cloud. The ap- 
parent stellar density, which makes a cluster detectable, results 
from both the true cluster density and the extinction variation 
toward its direction. Although a classical clustering analysis is 
biased in this configuration, it is possible to extract relevant 
information on clusters with a statistical analysis which takes 
the extinction into account (Carpenter, 2000; Cambresy et al., 
2002). 

2. Cluster mapping 

The Perseus molecular cloud has an elongated and filamentary 
structure, as shown in the extinction map presented in Fig. 1. 
This map is obtained by measuring the 2MASS H - K s color 
excess in cells with adaptive size. The cell size is chosen to 
contain only 3 stars, and the full map is then convolved by an 
adaptive kernel to obtain a final uniform spatial resolution of 
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Fig. 1. Perseus extinction map from 2MASS H - K s color excess. The 
spatial resolution is 4' and isocontours are for A v =2, 4, 8, 16 mag. The 
circle position and radius are set to contain all the confirmed members 
ofIC348. 



4'. This technique reduces the non-linear error because of the 
extinction variations at a smaller scale than the map resolution 
(Cambresy et al., 2005). The known members of IC 348 were 
removed before the extinction map was built. 

The final map is very similar to the 13 CO map (Padoan 
et al., 1999) and to the extinction map presented on the 
COMPLETE webpage 1 with the NICER method (Lombardi & 
Alves, 2001). An extinction estimate from reddening or from 
star counts relies on the assumption that all stars are behind the 
obscuring cloud and that the stellar population is homogeneous 
over the whole field. When a young stellar cluster is embed- 
ded in a molecular cloud, none of the conditions is fulfilled. 
The color and the star count maps react differently to the clus- 
ter contamination. Cambresy et al. (2002) proposed using this 
property to identify clusters by comparing both maps. The next 
3 sections describe the technique to build the color map and the 
count map and to obtain the cluster morphology. 

In the following, we focus on a 2° x 2° field centered on IC 
348atglon= 160.38°, glat= -17.78°. 

2.1. Reddening 

The technique is based on the pioneering work by Lada et al. 
(1994), who propose a statistical method to build extinction 
maps from near-infrared color excesses. The specificity of our 
method is to use adaptive cells that contain a fixed number of 
stars, in order to increase the spatial resolution and to limit the 
bias due to the inhomogeneous distribution of stars that tend to 



concentrate in the less obscured part of a cell. A detailed de- 
scription of this bias can be found in Cambresy et al. (2005). 
The visual extinction was obtained through the 2MASS H - K s 
color excess as follows: 
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where the indices obs and int stand for observed and intrinsic, 
respectively. Here, (H - K s )\ nl = 0.12 + 0.1 mag from the 
region outside the Perseus cloud. We used the extinction law 
from Rieke & Lebofsky (1985), for which A H /A V = 0.175 and 
AkJAv = 0.112, yielding the color excess to visual extinc- 
tion conversion formula A v = 15.87 x E H - Ks . The presence of 
a young embedded cluster will affect the extinction, because 
1) YSOs are embedded in the cloud and suffer only a partial 
obscuration, and 2) YSOs may have different colors than field 
stars, with an excess at K s due to an accretion disk for those 
objects that are younger than about 2 Myr. In the former case, 
the extinction will be underestimated, in the latter it will be 
overestimated. The effect does not usually stand out when ex- 
amining the reddening map because no anomaly (e.g. strong 
gradient or discontinuity) appears there. Another independent 
extinction estimation is required to detect the cluster contami- 
nation. 

2.2. Star counts 

In the star count method, the decrease in the apparent stellar 
density is interpreted as a foreground obscuration. It is the his- 
torical method for mapping the interstellar extinction (Wolf, 
1923; Bok, 1956). As mentioned above, we used an adaptive 
technique with the exact same cells as for the color excess map- 
ping. The visual extinction is obtained by assuming the appar- 
ent variations of the stellar density are due to foreground ob- 
scuration: 
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where a is the slope of the K s luminosity function (KLF), D 
the star density in each cell, and D re f the density in a reference 
field. Although the reference star density formally depends on 
the position (mainly on the galactic latitude), a constant value 
is correct for small areas of about 1 deg 2 . When a cluster con- 
taminates the field, the stellar density is locally enhanced, and 
the extinction is biased toward low values. If the local density 
is higher than the average density over the field (D re f), which is 
true if the extinction is small, it gives a negative number for the 
extinction. Obviously, Eq. 2 is no longer correct, because the 
measured density actually includes the background stars and 
the cluster stars : D = D\, + D c \ 

2.3. Comparison of the two maps 

Comparison of the maps requires the two calibrations to be 
consistent. The spatial resolution is identical in both maps by 
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Fig. 2. K s luminosity function of the 297 stars of the cluster with avail- 
able K s magnitude (Luhman et al., 2003, 2005a). 

construction. The relative calibration of the two maps is ob- 
tained using an arbitrary value for D re f in Eq. 2 and by forcing 
the median of AA V = A™ 1 - A™ 1 to zero. The absolute calibra- 
tion of the extinction maps is not relevant for the cluster study. 

The difference AA V should be a structureless map, if the 
field is free of any cluster contamination. Otherwise, the cluster 
morphology will be revealed by the structures contained in the 
AA V map because D c \ > 0, meaning A c ° l > A™'. 



3. Application to the IC 348 young cluster 

3.1. Data 

The data were extracted from the 2MASS All-Sky Point Source 
Catalog (Cutri et al., 2003) for the region defined by 158.6° < 
glon < 162.2° and -19.2° < glat < -16.3°. The list of the 304 
confirmed members of IC 348 was provided by Luhman et al. 
(2003) and Luhman et al. (2005a), 278 of which are present 
in 2MASS. Figure 2 shows the K s luminosity function of the 
full member list. The distribution peaks at K s « 1 1 .5 mag with 
more than 75% of the sources brighter than K s = 13 mag. On 
the other hand, faint stars dominate the background star pop- 
ulation with a number of sources proportional to \0 axK \ with 
a = 0.31. While a magnitude cut at K s < 13 mag keeps 75% of 
the cluster members, it removes 75% of the background stars. 
In order to increase the signal-to-noise ratio in the AAv map, 
we applied this cut before building the reddening and star count 
maps. 
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Fig. 3. Surface density map of IC 348 obtained from stars brighter 
than K s = 13 mag. The map shows the region 2<x = 7.7 mag above the 
noise. 
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Fig. 4. Same as Fig. 3 but derived from a catalog filtered from all 
known IC 348 members. These members are overlaid on the map as 
crosses. Circles at 5', 10', and 25' radius define the core (Lada & Lada, 
1995), the halo (Muench et al., 2003), and the cluster boundary sug- 
gested from our study. 



3.2. Results 

Following the method of Sect. 2, the extinction was from adap- 
tive cells that contain 3 stars. The reddening and the count maps 
were then smoothed with a Gaussian kernel to reach the final 
uniform resolution of 8'. We stress that it requires a convolu- 
tion by an adaptive kernel with FWHM = (8 2 - FWHM 2 rig ) 1/2 , 
where FWHM or ; g is the resolution in arcmin of the original 
maps. The difference between the maps from Eqs. 1 and 2 is 



presented in Fig. 3. As expected a peak appears at the IC 348 
position. 

The very same process was repeated on the 2MASS cat- 
alog filtered from all known IC 348 members. The resulting 
AAy map is presented in Fig. 4 with the cluster members 
overlaid as crosses. Significant structures remain on the map, 
which strongly supports the idea that the cluster extends be- 
yond the halo defined by Muench et al. (2003) at 10.33' and 
beyond the ROSAT X-ray detection boundary at 15' (Preibisch 
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et al., 1996). The map exhibits a faint ring around the cen- 
tral part of the cluster and several substructures. The farthest 
clump, located at / = 160.26 b = -18.42, is 38' south of 
the IC 348 center and seems linked to the main cluster by 
a ridge. It contains an infrared source reported by IRAS and 
MSX, MSX6C G160.2784-18.4216 (Kraemer et al., 2003). It 
is, however, likely that this structure corresponds to another 
star-formation site in Perseus. Further analysis needs to be car- 
ried out to conclude whether it is related to IC 348. Another 
clump located to the north-west of IC 348 has the particular- 
ity of being within the halo radius (Muench et al., 2003) of 
IC 348, which means this region has already been studied, al- 
though several bright members were not identified. 

Figures 3 and 4 were obtained by subtracting two extinction 
maps, making the pixel unit to be magnitude. It has no physical 
meaning in this context, except that it goes with the star den- 
sity. A straightforward way to convert these magnitudes into 
star density would require the real extinction to be known, but 
both the reddening and the count maps are contaminated by the 
YSOs. Fortunately, the outer region of IC 348 has a low star 
density, since it is the first time the excess has been detected. 
It means the sample must be dominated by background stars, 
which makes the H - K s color excess an acceptable estima- 
tion of the extinction. We stress that the median color is used 
within cells rather than the mean, which ensures a low sensi- 
tivity to the YSOs contamination as long as it is less than 50%. 
Assuming the color-based extinction map is correct, Fig. 4 was 
converted in surface density as follows: 

D c i = cx(l0 flAA *' - l) (3) 

where c is obtained using the knowledge that the difference 
between N c \ in Figs. 3 and 4 must be exactly the number 
of confirmed members that were removed from the original 
2MASS catalog to make the second figure (i.e. 212 stars for 
K s < 13 mag). We finally conclude that 40 + 6 stars with 
K s < 13 mag are responsible for the excess seen in Fig. 4. 
The quoted uncertainty is the statistical uncertainty associated 
with the Poisson law followed by the star distribution. A total 
number of 350 stars brighter than 13 mag at K s was found in 
the area, meaning an 11% level contamination by the cluster, 
which justifies our choice of relying on the color median as an 
extinction estimator. Following Eq. 3, the southern clump cen- 
tered on MSX6C G160.2784-18.4216 would contain about 10 
YSOs brighter than K s < 13 mag. 

The color-color diagram (Fig. 5) for the stars located in the 
region defined in Fig. 4 shows almost no star with a H - K s 
excess. This is also true for the confirmed IC 348 members, 
which indicates that most YSOs have already lost their inner 
accretion disk that is responsible for the infrared excess at 2 
fim. Individual identification is therefore not possible with the 
2MASS photometry. 

4. Conclusions 

The statistical comparison of the stellar color and density per- 
mit building the surface density map of embedded clusters. We 
applied our technique to IC 348, using two different sample cat- 
alogs, one with all 2MASS sources and the other one filtered 
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Fig. 5. Color-color diagram for stars with K s < 13 mag in the area 
defined by Fig. 4 (diamonds) and for confirmed IC 348 members 
(crosses). The dashed lines are parallel to the reddening vector. 

from all known associated YSOs. As expected, the first map 
shows a peak at the cluster location, but it is more surprising 
to find evidence of structures in the filtered map, too. These 
structures are distributed around IC 348, suggesting they are 
actually an extension of about 40 new members for the well- 
known cluster. This result encourages more detailed study of 
the structure of IC 348 and, in particular, of its dynamical evo- 
lution, as it created these clumps in the cluster. As no infrared 
excess was found in 2MASS, we cannot identify the potential 
new members individually. Longer wavelength follow-up with 
Spitzer may exhibit the presence of a residual disk, but X-ray 
observations with XMM-Newton are probably more appropri- 
ate to finding the older pre-main-sequence stars in the outer re- 
gion of IC 348. Spectroscopic follow-up remains the most suit- 
able way to obtain the final diagnostic for membership. Proper 
motion measurements could be an indicator of the membership 
and would be of great interest for the kinematic study of the 
cluster. 
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